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ABSTRACT 

In analogy to f(R) theory, recently a new modified gravity theory, namely the so-called f(T) 
theory, has been proposed to drive the current accelerated expansion without invoking dark energy. 
In the present work, by extending Bisabr's idea, we try to constrain f(T) theories with the varying 
fine structure "constant", a = e 2 /he. We find that the constraints on f(T) theories from the 
observational Aa/a data are very severe. In fact, they make f(T) theories almost indistinguishable 
^ ! from ACDM model. 
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I. INTRODUCTION 

The current accelerated expansion of our universe [l[ has been one of the most active fields in modern 
cosmology since its discovery in 1998. This mysterious phenomenon could be due to an unknown energy 
component (dark energy) or a modification to general relativity (modified gravity) [H-Q • The well-known 
modified gravity theories are, for examples, f(R) theory, scalar-tensor theory (including Brans-Dicke 
theory), braneworld scenarios (such as DGP, RSI and RSII), f(G) theory (G is the Gauss-Bonett term), 
Horava-Lifshitz theory, MOND and TeVeS theories. We refer to e.g. PHD, [5(| for some reviews. 

Recently, a new modified gravity theory, namely the so-called f(T) theory, attracted much attention 
in the community, where T is the torsion scalar. It is a generalized version of the so-called teleparallcl 
gravity originally proposed by Einstein [H, Hj] . In teleparallel gravity, the Weitzenbock connection is used, 
rather than the Levi-Civita connection which is used in general relativity. Following 0, H|, here we 
briefly review the key ingredients of teleparallel gravity and f(T) theory. We consider a spatially flat 
Friedmann-Robertson- Walker (FRW) universe whose spacetime is described by 

ds 2 = -dt 2 + a 2 (t)dx 2 , (1) 

where a is the scale factor. The orthonormal tetrad components e^x 1 *) relate to the metric through 

g^u = %-e^4 , (2) 

where Latin i, j are indices running over 0, 1, 2, 3 for the tangent space of the manifold, and Greek 
H, v are the coordinate indices on the manifold, also running over 0, 1, 2, 3. In teleparallel gravity, the 
gravitational action is 

*r = i/^HT, (3) 



2k 2 j 

where k 2 = 8irG, and |e| = det (ejj = y/—g ■ The torsion scalar T is given by 

r = (4) 

where 

TV = ~4(d»<-d^ , (5) 



K^ p = -- (T^p - T v % - T/") , (6) 



Sr = \ + 8p* v e ~ 5 V pT e » e ) ■ (7) 

For a spatially flat FRW universe, from Eqs. and ([TJ, one has 

T = —QH 2 , (8) 

where H = a /a is the Hubble parameter (a dot denotes the derivative with respect to cosmic time t). So, 
one can use T and H interchangeably. In analogy to f(R) theory, one can replace T in the gravitational 
action © by a function f(T) (see however Q). In f(T) theory, the modified Friedmann equation and 
Raychaudhuri equation are given by @, ||| 

12H 2 f T + f = l6TrGp, (9) 

48H 2 f TT H - fx (l2H 2 + m) - f = 16nGp , (10) 



where /t = df /dT, and p, p are the total energy density and pressure, respectively. In an universe with 
only dust matter, p = p rn = and p = p rn . From Eqs. @ and (JTDJ), one can find that the effective dark 
energy density and pressure from torsion are given by [3, HI, flO| 

Pde = (6H 2 - f - \2H 2 f T ) , (11) 

Pde = -Pde ~ ^ (l2H 2 f TT -f T + l)H. (12) 

Obviously, if f(T) = T + const., f(T) theory reduces to the well-known ACDM model. 
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In fact, /(T) theory was firstly used to drive inflation by Ferraro and Fiorini [THOU]- Later, Bengochea 
and Ferraro Q, as well as Linder Q, proposed to use f(T) theory to drive the current accelerated 
expansion without invoking dark energy. Soon, many works followed. For examples, Myrzakulov [l3T | and 
Yang [l4[ proposed some new f(T) forms; Bengochea [l5j], Wu and Yu [l6| considered the cosmological 
constraints on f(T) theories by using the latest observational data; Wu and Yu [13] also considered the 
dynamical behavior of f(T) theory; Dent et al. [l8j], Zheng and Huang [Ti| considered the cosmological 
perturbations and growth factor in f(T) theories; Wu and Yu [2(|, Bamba and Geng [2l| discussed 
the equation-of-state parameter (EoS) crossing the phantom divide in f(T) theories; Zhang et al. [13] 
discussed the dynamical analysis of f(T) theories; Li, Sotiriou and Barrow [H| considered the large- 
scale structure and local Lorentz invariance in f(T) theory; Deliduman and Yapiskan [24| discussed the 
relativistic neutron star in f(T) theory; Cai et al. [25j | considered the matter bounce in f(T) theory; 
Wang 26] discussed the static solutions with spherical symmetry in f(T) theories. We further refer to 
e -g- [23] for some relevant works. 

In the literature, the observational constraints on f(T) theories [HI, EH were obtained mainly by using 
the cosmological data, such as type la supernovae (SNIa), baryon acoustic oscillation (BAO), and cosmic 
microwave background (CMB). In the present work, we instead try to constrain f(T) theories with the 
varying fine structure "constant", a = e 2 /he. In Sec. UH we briefly review the observational constraints 
on the temporal variation of the fine structure "constant" a. In Sec. IIII1 we briefly introduce the idea to 
constrain f(R) theories with the varying a, which was proposed by Bisabr [28]. In Sec. IIV1 we extend 
Bisabr's idea to f(T) theories, and consider the corresponding constraints from the temporal variation 
of the fine structure "constant". Finally, some brief concluding remarks are given in Sec. fVl 



Aa/a 


Redshift 


Observation 


Ref. 


<io~ 2 


10 10 - 10 s 


BBN 


[33, 34] 


< 1(T 2 


10 3 


CMB 


[34] 


< 1(T 6 


3-0.4 


quasars 


[30, 31 , 35 , 36} 




0.45 


meteorite 


[37] 


<io~ 7 


0.14 


Oklo 


[38] 



TABLE I: The observational constraints on Aa/a. 



II. OBSERVATIONAL CONSTRAINTS ON THE TEMPORAL VARIATION OF 
THE FINE STRUCTURE "CONSTANT" 

Motivated by the well-known large number hypothesis of Dirac proposed in 1937 29], the varying 
fundamental "constants" remain as one of the unfading subjects for decades. Among the fundamental 
"constants", the most observationally sensitive one is the electromagnetic fine structure "constant", 
a = e 2 /he. Since about 12 years ago, this subject attracted many attentions, mainly due to the first 
observational evidence from the quasar absorption spectra that the fine structure "constant" might change 
with cosmological time [13, Hlj] . 

Subsequently, many authors obtained various observational constraints on the temporal variation of the 
fine structure "constant" a. In the literature, it is convenient to introduce a quantity Aa/a = (a— ao)/ao, 
where the subscript "0" indicates the present value of the corresponding quantity. Obviously, Aa/a is 
time-dependent. A brief summary of the observational constraints on Aa/a can be found in e.g. [32l ]. The 
most ancient constraint comes from the Big Bang Nucleosynthesis (BBN) [33|,[34j], namely, | Aa/a\ < 10~ 2 , 
in the redshift range z = 10 10 — 10 8 . The next constraint comes from the power spectrum of anisotropy in 
the cosmic microwave background (CMB) (34|, i.e., (Aa/a] < 10~ 2 , for redshift z ~ 10 3 . In the medium 
redshift range, the constraint comes from the absorption spectra of distant quasars [13, HH, HI, H3|- 
Since the results in the literature are controversial, it is better to consider the conservative constraint 
|Aa/a| < 10~ 6 [llljjn the redshift range z = 3 — 0.4. From the radioactive life-time of 187 Re derived from 
meteoritic studies [33, the constraint is given by |Aa/a| < 1CP 7 for redshift z = 0.45. Finally, from the 
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Oklo natural nuclear reactor (38[, it is found that \Aa/a\ < 10 7 for redshift z — 0.14. For convenience, 
we summarize the above constraints in Table [H which will be used in the followings. 



III. THE IDEA TO CONSTRAIN f(R) THEORIES WITH VARYING ALPHA 
A. Varying alpha driven by a general scalar field 



Noting that a = e 2 /he, a varying a might be due to a varying speed of light c [39l44ll ]. while Lorentz 
invariance is broken. The other possibility for a varying a is due to a varying electron charge e. In 1982, 
Bekenstein proposed such a varying a model (42j, which preserves local gauge and Lorentz invariance, and 
is generally covariant. This model has been revived a nd g eneralized after the first observational evidence 
of varying a from the quasar absorption spectra [3(| l3ll | This is a dilaton theory with coupling to the 
electromagnetic F 2 part of the Lagrangian, but not to the other gauge fields. Later, the Bekenstein-type 
varying a model has been generalized by replacing the dilaton with a cosmological scalar field. Further, 
the coupling between the scalar field and the electromagnetic field could also be generalized. In fact, 
the var ying a models driven by quintessence have been extensively investigated in the literature (see 
e.g. [32, 43 4GJ). The varying a driven by phantom has been considered in the BSBM model [47| while 
its model parameter u> is negative. The special case of varying a driven by fc-essence whose Lagrangian 
£(X, 4>) = X n - V{4>) has been considered in e.g. 0. The varying a driven by Dirac-Born-Infeld scalar 
field has also been discussed in [isj . 

Following [H, |H, |4|| , the relevant action in Einstein frame is generally given by 

S = S g + J dtxy/^gCt-^j d 4 x^9B F (^)F^F^ +S m , (13) 

where S g is the gravitational action in Einstein frame; F^ are the components of the electromagnetic 
field tensor; S m is the action of other matters:/^ is the Lagrangian of the scalar field (j>. Noting that 
Bp takes the place of e~ 2 in Eq. (|13l) actually @, H^], one can easily see that the effective fine structure 
"constant" a = e 2 /hc is given bv |32tl43l] 



«o 



Thus, we find that 



B F ((f>( X ,t)) 
Aa a — aQ 1 



a ao Bf{4>) 



(14) 



(15) 



It is worth noting that the present value (at redshift z = 0) of the coupling Bp should be 1 by definition. 
If Bp(z = 0) = Bpo 7^ 1, we can normalize it through rescaling the electromagnetic field, namely 



Fftu -> VBpoF^ , while Bp -> — — . (16) 

-DF0 



Thus, we have 

Aa Bpo 
a Bp 



1. (17) 



B. Varying alpha in f(R) theories 

Following Bisabr's idea (28|, here we briefly show why the fine structure "constant" should be varying 
in f(R) theories. As is well known, in Jordan frame the action of f(R) theories with matters (including 
electromagnetic fields here) reads (see e.g. 0-Q) 



S =2 



^ J d 4 x^f(R) + J d 4 x C M (g^ M ), (18) 
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where R is the Ricci scalar, and Cm is the matter Lagrangian depending on g^ u and matter fields ^ m 
(including electromagnetic fields here). It is well known that f{R) theory can be equivalent to scalar- 
tensor theory 043] • Applying the conformal transformation 

Of 

9^ = V , n 2 = T = f R = ^ , (19) 
and introducing a new scalar field <j) defined by 

nip = \/s/2 In J- , (20) 
one can rewrite the action (fT8| to the one in Einstein frame namely 

R 1 



S = Id xy'—g 



+ / d i xC M {^ 1 ^)~9^^M) , (21) 



where a tilde denotes quantities in Einstein frame, and the potential of scalar field is given by 

yw-^f- (22) 

From Eq. (|2ip , it is easy to see that the matter fields (including electromagnetic fields here) are inevitably 
coupled with the scalar field <fi in Einstein frame. So, a should be varying. Noting that y/—g = fi~ 4 \/— ~g, 
one can easily find that the coupling in Eq. fj X3[) is given by [11] 

B F =T- 2 = f R 2 . (23) 

From Eq. (|T5|) . the variation of the fine structure "constant" can be described by [28[ 

a 

As mentioned above, if Jr{z = 0) = fno =/= 1, we can normalize it through Eq. (|16[) , and then 
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1. (25) 



Aa = (h_ 

a V /ho , 

In [2^, Bisabr discussed the observational constraints on f(R) theories with the temporal variation of 
the fine structure "constant", and found that the corresponding constraints are fairly tight. We refer to 
the original paper [28] for details. 

IV. f(T) THEORIES AND VARYING ALPHA 

As mentioned in Sec.|H f(T) theories are proposed in analogy to f(R) theories. So, we extend Bisabr's 
idea [28[ to constrain f(T) theories also with the temporal variation of the fine structure "constant". 

A. Varying alpha in a general f(T) theory 

As mentioned in [1, H3], f(T) theory can also be equivalent to scalar-tensor (torsion) theory. In Jordan 
frame, the relevant action with matters (including electromagnetic fields here) reads 

S=^Jd i x \e\ f(T) + J d 4 xC M (4, *m) . (26) 

Similarly, applying the conformal transformation 

V = "V o e*=fie* n 2 =T = f T , (27) 



G 



one can also rewrite the action ([26]) to the one in Einstein frame with a new scalar field <f> [1, E3 ■ Similar 
to the case of f(R) theories, the matter action in Einstein frame is given by 



§ M = J d 4 x£ M (F- 1/2 (4)el,*M 



(28) 



Again, the matter fields (including electromagnetic fields here) are inevitably coupled with the scalar field 
(ft in Einstein frame. So, a should be varying. Noting that |e| = f2 _4 |e|, one can similarly find that 



Bp — T frp ■ 

From Eq. (P~5]) . the variation of the fine structure "constant" can be described by 



Aa 
a 



— = It ~ 1 



As mentioned above, if fr{z = 0) = fra ^ 1, we can normalize it through Eq. (|16[) , and then 

2 



Aa 

a 



fro 



1 . 



(29) 



(30) 



(31) 



In the followings, we will consider two concrete f(T) theories, namely, f(T) = T + /i(— T) n and 
f(T) = T — \xT (l — e^ T °/ T ) which are the most popular f(T) theories discussed extensively in the 
literature (see e.g. 0, H EE El ) . 




O 




FIG. 1: log \Aa/a\ as a function of redshift z for f(T) = T + u(-T) n with n m0 = 0.272 and n = 0.04 (solid curve), 
f2 m o = 0.272 and n — 0.26 (thick long-dashed curve), Q m o = 0.272 and n — —0.29 (thick short-dashed curve), 
f2 m o = 0.240 and n = 0.04 (thin long-dashed curve), Q. m o = 0.308 and n = 0.04 (thin short-dashed curve). Right 
panel is an enlarged part of left panel. Only the curves not overlapping the gray areas are phenomenologically 
viable. See text for details. 



B. f(T) =T + At(-T)' 1 

At first, we consider the case of f(T) = T + /i(— T) n , where \x and n are both constants. This is the 
simplest model, and has been considered in most papers on f(T) theory. Obviously, if n = 0, it reduces 
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FIG. 2: log|Aa/a| as a function of redshift z for f(T) = T+n(-T) n with a fixed Q m0 = 0.272, and n = ±f.8xl(r 7 
(solid curve), n = ±10" 6 (long-dashed curve), n — ±1CT 8 (short-dashed curve). Only the curves not overlapping 
the gray areas are phenomenologically viable. See text for details. 



to ACDM model. Substituting it into the modified Friedmann equation one can easily find that /i is 
not an independent model parameter, namely @, fisj ] 



1 f^mO ( nT T2\ 1 - n 1 — ^mO / m \l—n 

— { 6H o) = — T (-To) 



2n 



2n-l 



where Q TO rj = 8ttG p m o / (3Hq) is the present fractional energy density of dust matter. So, we have 

2n — 1 V Jq 



and then 



/t = 1 



n(l - Omo) c,2(n-l) 



1 - 271 



/to — 1 



n(l - Qroo) 
1 - 2n 



(32) 



(33) 



(34) 



where E 12 = T/T = H 2 /H%. Substituting Eq. flMJ) into Eq. one can finally obtain the explicit 

expression of Aa/a. In order to compare it with the observational data, we need to know E(z) as a 
function of redshift z. Substituting f(T) = T + fj,(—T) n and Eq. (|3"2"|) into the modified Friedmann 
equation ©, we find that 



e 2 = n m0 {i + z) 3 + (i-n m0 )E 2 



(35) 



Obviously, if n — 0, it reduces to the one of ACDM model. If il m o and n are given, we can numerically 
solve Eq. (135]) and obtain E 2 {z) as a function of redshift z. Thus, Aa/a is on hand. 

Next, we compare Aa/a with the observational data. In fact, as shown in e.g. [32l.l48|. the constraints 
from the first two rows (at very high redshift) in Table U are very weak. Therefore, we only consider the 
last three rows (at low redshift) in Table U (and hence the radiation can be safely ignored). Note that 
in [lfij], this f(T) = T + n(—T) n model has been constrained by using the latest cosmological data, i.e., 
557 Union2 SNIa dataset, BAO, and shift parameter from WMAP7. The corresponding 2er results are 
given by [l6| 



(36) 
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At first, we try to see whether Aa/a with the best-fit parameters of [1 61 ] and the corresponding 2a edge can 
simultaneously satisfy the observational constraints in TableQ] In Fig. [I] we plot log |Aa/a| as a function 
of redshift z for f(T) = T + ju(-T)™ with Q m0 = 0.272 and n = 0.04 (solid curve), £l m0 = 0.272 and 
n = 0.26 (thick long-dashed curve), fl m o — 0.272 and n — —0.29 (thick short-dashed curve), fl m o — 0.240 
and n = 0.04 (thin long-dashed curve), f2 m o = 0.308 and n = 0.04 (thin short-dashed curve), where log 
indicates the logarithm to base 10. Obviously, one can see that Aa/a with the best-fit parameters of [l|| 
and the corresponding 2a edge cannot satisfy the observational constraints in Table [TJ In addition, from 
Fig. [1] we find that the influence from n to Aa/a is significantly larger than the one from f2 m o- Thus, 
fixing O to q = 0.272, we try various n to find in which cases all the observational constraints in Table U 
could be simultaneously satisfied. From Fig. [3J it is easy to see that they can be all respected only for 

\n\ < 1.8 x 10~ 7 . (37) 

This is the constraint on f(T) = T + fi(—T) n theory from the observational Aa/a data. It is a very 
severe constraint in fact. Noting that f(T) = T + /i(—T) n — » T + const, when n — > 0, this f(T) theory 
becomes almost indistinguishable from ACDM model. 
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^ -4 

< 

o -6 




f(T) = T - fj,T(l - e? 1 " 11 ) 



0.5 1 1.5 2 2.5 3 
z 



0.5 1 1.5 2 2.5 3 



FIG. 3: log [Aa/a | as a function of redshift z for f(T) = T — fj,T f 1 — e BT ° /T ) with Q m0 = 0.272 and /3 = -0.02 
(solid curve), Q m0 = 0.272 and f3 = 0.29 (thick long-dashed curve), f2 m0 = 0.272 and f3 = -0.22 (thick short- 
dashed curve), QmO = 0.238 and /3 = —0.02 (thin long-dashed curve), fl m o = 0.308 and /3 = —0.02 (thin 
short-dashed curve). Right panel is an enlarged part of left panel. Only the curves not overlapping the gray areas 
are phenomenologically viable. See text for details. 



C. /(T)=T-/xT(l 



Here, we consider the case of f(T) = T—fiT (l — e /9T °/ T ) , where [i and /3 are both constants. Obviously, 
f(T) -» T+fif3T Q = T + const. when /3 — > 0, it reduces to ACDM model. Substituting it into the modified 
Friedmann equation one can easily find that /i is not an independent model parameter [8], El, i.e., 



/' 



1 — QmO 



1 - (1 - 2/3) e 
On the other hand, it is easy to obtain 

/ T = l- M + M (l-/3/i? 2 ) e ^ 2 , 



/to = 1 - M + M (1 - P) 



(38) 



(39) 
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f(T) = T - jiT(l - e f3T o /T ) 
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FIG. 4: log |Aa/a| as a function of redshift z for f(T) = T — /j,T (l — e 8To/T j with a fixed Q m0 = 0.272, and 

P = ±2.3 x 1CT 7 (solid curve), ft = ±10 6 (long-dashed curve), /3 = ±10 7 (short-dashed curve). Only the curves 
not overlapping the gray areas are phenomenologically viable. See text for details. 



where E 2 = T/T = H 2 /H%. Substituting Eq. fl3H]) into Eq. dSU) , one can finally obtain the explicit 
expression of Aa/a. In order to compare it with the observational data, we need to know E(z) as a 
function of redshift z. Substituting /(T) = T—/iT (l — e /3T °/ T ) into the modified Friedmann equation (|9]), 



we find that 



E 2 = n m0 (l ± zf + fiE 2 



(40) 



If (3 — > 0, we have fip -> 1 — tt m0 from Eq. ((38)) . and hence Eq. (|40| reduces to the one of ACDM model 



z. 



If 57 m o and /3 are given, we can numerically solve Eq. (|4U| and obtain E 2 (z) as a function of redshift 
Thus, Aa/a is on hand. 

In [lfij], this f(T) = T — /iT (l — e /3T °^ T ) model has also been constrained by using the latest cosmo- 
logical data, i.e., 557 Union2 SNIa dataset, BAO, and shift parameter from WMAP7. The corresponding 
2cr results are given by 16] 

O m0 = 0.272+0;°^ , (3 = -0.02±^i . (41) 

Again, we try to see whether Aa/a with the best-fit parameters of [l6[ and the corresponding 2a edge 
can simultaneously satisfy the observational constraints in Table |U In Fig. [21 we plot log |Aa/a| as a 
function of redshift z for f(T) = T — (j,T (l — e /3T °/ T ) with Q m0 = 0.272 and (3 = -0.02 (solid curve), 
n m0 = 0.272 and (3 = 0.29 (thick long-dashed curve), ft m0 = 0.272 and /3 = -0.22 (thick short-dashed 
curve), fl m0 = 0.238 and (3 = -0.02 (thin long-dashed curve), Q m0 = 0.308 and /3 = -0.02 (thin 
short-dashed curve). Obviously, one can see that Aa/a with the best-fit parameters of [16( and the 
corresponding 2a edge cannot satisfy the observational constraints in Table HI In addition, from Fig. |3l 
we find that the influence from (3 to Aa/a is significantly larger than the one from J7 m o- Thus, fixing 
^mo = 0.272, we try various (3 to find in which cases all the observational constraints in Table U could be 
simultaneously satisfied. From Fig. HI it is easy to see that they can be all respected only for 

\fi\ < 2.3 x 10~ 7 . (42) 

This is the constraint on f(T) = T — /iT (l — e^ T °^ T ) theory from the observational Aa/a data. It is a 
very severe constraint in fact. Noting that f(T) = T — jiT (l — e /3T °/ T ) — > T + const, when j3 — > 0, this 
f{T) theory becomes almost indistinguishable from ACDM model. 
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V. CONCLUDING REMARKS 

In analogy to f(R) theory, recently f(T) theory has been proposed to drive the current accelerated 
expansion without invoking dark energy. In the literature, the observational constraints on f(T) theories 
were obtained mainly by using the cosmological data, such as type la supernovae (SNIa) , baryon acoustic 
oscillation (BAO), and cosmic microwave background (CMB). In the present work, by extending Bisabr's 
idea (28|, we instead try to constrain f(T) theories with the varying fine structure "constant", a = e 2 /he. 
We found that the constraints on f(T) theories from the observational Aa/a data are very severe. In 
fact, they make f(T) theories almost indistinguishable from ACDM model. 

Some remarks are in order. Firstly, in this work we only considered a spatially flat FRW universe. 
This is mainly motivated by the well-known inflation scenario and the very tight observational constraint 
on the spatial curvature term from WMAP7 data [5lj |. namely, ilk = — 0.0057^0 0068- Obviously, when 
a non- vanishing spatial curvature term is allowed, the observational constraints on f(T) theories from 
the varying fine structure "constant" should be relaxed, since the number of free model parameters 
is increased. However, we can expect that the situation of f(T) theories are still not improved, due 
to the very narrow range of the constrained Ofc. For instance, even the constraints on n or /3 could be 
greatly relaxed from 0(1O -7 ) to 0(1O -4 ) (say), the corresponding f(T) theories are still indistinguishable 
from ACDM model. Secondly, in fact the fine structure "constant" might be not only time-dependent 
but also space-dependent (see e.g. [52j|). Of course, the space-dependent fine structure "constant" is 
still in controversy. On the other hand, the space-dependent fine structure "constant" might invoke an 
inhomogeneous scalar field cj). As is shown in Sees. IIIIBI and IIV Al the f(R) and f(T) theories in a 
homogeneous and isotropic FRW universe could not lead to a space-dependent fine structure "constant" , 
and hence they are not constrained by the possibly spatial variation of the fine structure "constant" . 
So, in this work we only considered the time-dependent fine structure "constant" for simplicity. Thirdly, 
it is worth noting that f(T) gravity does not generally preserve the local Lorentz invariance and any 
theory of f(T) gravity is always built on a local frame which is chosen on a specific spacetime point. As 
consequences, it is not always applicable for the conformal transformation in f(R) gravity to be used 
in f(T) gravity, and it is quite unclear how to explicitly define an Einstein or a Jordan frame in f(T) 
gravity. Fortunately, we could use these conceptions in the homogenous and isotropic FRW background 
which is quite a special case (we thank the referee for pointing out this issue). Finally, we note that 
the two concrete f(T) theories considered in this work contain only a single free parameter, which are 
very simple cases. In fact, it is reminiscent of the case of f(R) theory. In the beginning, the forms of 
/(i?) are also very simple, such as the types of 1/R or R n [IHJ]. Later, these simple f(R) forms have 
been easily ruled out by the cosmological observations and the local gravity tests. After several years, 
the viable f(R) forms which can satisfy all the cosmological observations and the local gravity tests, e.g. 
Hu-Sawicki (53|, Starobinsky [5J] and Tsujikawa (55[, have been toughly earned [JH|- These three viable 
f(R) forms are all complicated and delicate, and they contain two or more free parameters. Similarly, we 
expect that the viable f(T) forms, which satisfy all the cosmological observations, the local gravity tests 
and the constraints from the varying fine structure "constant" , could be constructed with tough efforts 
in the future. Of course, it is easy to anticipate that they are also complicated and delicate, and contain 
many free parameters. In addition, at that time it is also interesting to see what kind of dark energy 
models could be mimicked by the viable f(T) theories. 
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